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Abstract 



Higher order derivative corrections to the Einstein-Maxwell action are considered and 
an explicit form is found for the corrections to the entropy of extremal black holes. We 
speculate on the properties of these corrections from the point of view of small black 
holes and in the case when the classical black hole potential exhibits flat directions. A 
particular attention is paid to the issue of stability of several solutions, including large 
and small black holes by using properties of the Hessian matrix of the effective black hole 
potential. This is done by using a model independent expression for such matrix derived 
within the entropy function formalism. 



1 Introduction 



Black hole physics provides us with a vast variety of phenomena for testing underlying ideas of 
theoretical physics. This explains the constant attention to this topical issue, as well as the burst 
of interest any time a new concept emerges in this area of research. Superstring/M-theory has 
enriched physics with new ideas and currently is the main subject of interest in mathematical 
physics. It allowed to intertwine different areas of theoretical physics, e.g. gravity and conformal 
field theories. 

There is a well known correspondence between black hole mechanics and thermodynamics [1] 
that relates geometrical characteristics of a black hole to thermodynamical ones. One of them, 
the subject of our interest here, - entropy - can be calculated using the Wald formula [2]. To 
be able to interpret this quantity as a genuine entropy, it should be confirmed by statistical 
physics calculations. This was accomplished by performing counting of the microstates [3i |4] in 
the "classical" limit. At this point, there arises a question of comparison between microstate 
and "macrostate" entropies beyond the classical approximation. 

For extremal black holes the method of calculation of macrostate entropy was proposed 
in [5]. It relies on the classical Einstein-Maxwell action and reduces the problem to finding 
a black hole potential which encodes all the information needed to obtain the entropy, which 
is given by the value of the black hole potential at its critical point. This method allows one 
to trace back the evolution of the scalar fields in the whole space, but might become difficult 
when considering higher order corrections which originate from the string coupling and the a ' 
expansion of the superstring action. 

For this purpose another approach was proposed in J6j. It deals with horizon values of the 
fields present in the theory by means of the so-called entropy function. 

Taking into account terms coming from the superstring action (such as the Chern-Simons 
term) allows one to achieve a matching between black holes microstate and "macrostate" en- 
tropies [7]. 

Higher order corrections might be of interest from the classical point of view due to the 
following problems: 

1. small black hole - black hole with vanishing classical entropy and non- vanishing values 
of mass. Once classically one obtains a zero value of the entropy, one naturally poses the 
question what happens when taking into consideration higher order corrections as well. 

2. flat directions of the black hole potential and stability of the critical points. The presence 
of fiat directions reflects the symmetry of the scalar manifold. Therefore it is interesting 
to know how the symmetry gets modified in the presence of higher order corrections and, 
as a consequence, how fiat directions get distorted and whether the resulting critical point 
is stable or not. 

The problem of small black holes is studied from different points of view [5]. 

In order to investigate the influence of higher order corrections on the solutions to the 
attractor mechanism equations [j9j, we make use of the black hole potential approach and 
the entropy function formalism. We established a direct relation between these methods and 
found how the black hole potential is related to the entropy function. In addition, we derive a 
formula flTUl) for the Hessian matrix of the black hole potential completely within the entropy 
function formalism. A distinctive feature of this formula is that it can be easily used even 
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if deriving the expression of the black hole potential from the entropy function proves to be 
problematic. 

Previously, in the entropy function formalism not much attention was paid to the stability 
issues, due to the fact that a truncation over the "problematic" sectors (those where the fiat 
directions reside, e.g. the axionic one in the stu model) was always considered. Having at 
disposal the formula for the Hessian matrix allows us to study the stability of the complete, 
non truncated version of the theory. 

In the present work, we consider the most general case of higher order derivative corrections 
to the Einstein-Maxwell action involving both Riemann curvature tensor Rnupa and electro- 
magnetic field strength tensor F^^^ in any Lorentz covariant combinations. From dimensional 
analysis it follows that the ra*^ order correction should have the form {R^,^p^)"^{Fl^j^Y^''^~"^'^^^ 
with m = 0, . . . , n + 1. The indices are supposed to be contracted in all possible ways by means 
of a metric g^y. The number of such terms grows very fast with respect to n, but in the entropy 
function formalism they all lead to a much smaller number of independent combinations. In 
this way we succeed in finding the corrections to the black hole entropy induced by the higher 
order derivative corrections to action. 

As it is known, usually, in order to calculate up to n**^ order the value of a function at a 
critical point, one should know the solution to the criticality condition up to (n — 1)*'^ order. 
Examining the form of the corrections to the entropy, one might easily observe that there 
appears a problem starting from the first order. Namely, the first order correction to the 
entropy is supposed to be determined by the classical solution to the attractor mechanism 
equations. Once the classical black hole potential possesses fiat directions, not all moduli are 
determined by the attractor mechanism equations, though this ambiguity does not affect the 
classical value of the entropy. But in the first order this ambiguity shows up: the entropy 
becomes dependent on the undetermined moduli. Evidently this fact contradicts the second 
law of black hole thermodynamics and the attractor mechanism paradigm. Nevertheless, a 
thorough analysis of the perturbative corrections shows that there exists a sort of feedback 
from the first order solutions on the classical ones. In general, in order to calculate the entropy 
up to n^^ order, one should have the solution up to n*^ order. As it was mentioned above, n^^ 
order solution does not show up in the expression for the entropy (at least to n*^ order), but 
it might not exist in the presence of fiat directions. This is related to the fact that, in order to 
find the n**^ order solution, one has to solve a system of linear non-homogeneous equations with 
vanishing determinant. This system contains as parameters moduli which are not defined in 
the previous orders. Some of these parameters get fixed when requiring the system of equations 
to be consistent. In this way the equations to n*'^ order fix the solutions in previous orders. 
Exactly for this reason the attractor mechanism remains valid in the presence of higher order 
derivative corrections, as well. 

We considered in detail an example of a stu dyonic black hole in heterotic string theory 
compactified on or K3 x T^. In addition to previously known solutions [TldU], we found two 
new solutions. They are both non-BPS and one of them, moreover, corresponds to a state with 
vanishing central charge Z |TT]. In the classical limit the latter solution turns out to be stable. 
Therefore we dwell mainly on the stability of the non-BPS Z ^ solution, which classically 
has two vanishing and four positive eigenvalues. We find that, when the corrections are turned 
on, one of the previously vanishing eigenvalues remains zero, while the other becomes positive. 

We studied as well the small black hole limit and found that the two new solutions mentioned 
above do not allow for such a limit, while the previously known ones do. Therefore we studied 
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the stability of these solutions, as well, and found that for both of them in the small black hole 
limit the Hessian matrix has always at least one negative eigenvalue. 

The paper is organized as follows. In section [2] we establish a relation between the two 
methods for calculating the black hole entropy (j6]). The main result of this section is the 
formula (ITU]) yielding the Hessian matrix expressed in terms of the entropy function. In section [3] 
we introduce the higher order corrections to the action, discuss their explicit form and derive 
the expressions for the black hole effective potential f|T7j) and entropy fl26l) . up to the second 
order. We speculate on properties of the corrections to the entropy both in the case of small 
black holes and when the classical black hole potential has flat directions. In section H] we 
revisit a well known example of heterotic string theory compactified on or K3 x without 
neglecting axions. We demonstrate here that the corrections to the entropy ( l39l) do not depend 
on values of the scalars at infinity. We derive two new non-BPS solutions and investigate their 
stability, as well as the stability of the previously known ones. A special attention is paid to 
small black holes and the issue of their stability. 



2 Basics 



Since the paper [5] , it has been known how to calculate the entropy of a wide class of extremal 
black holes. There an Einstein-Maxwell model coupled to scalar fields was considered 
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(1) 



and it was shown that the entropy of a black hole in a static, spherically symmetric and 
asymptotically fiat space-time is completely determined by the so-called black hole potential, 
which for the model ([T]) has the form 



Vbh(0) 



-2^P^^ (1a) 
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The entropy within this framework is given by the value of the black hole potential at the 
horizon 

S = 7tVbhM, (3) 
where the moduli take critical values obtained from 
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This approach naturally allows us to examine the stability (which guarantees that the expression 
in ([3]) refers to the entropy of a physical object) of the critical points by checking the positive 
definiteness of the Hessian matrix 
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and it gives us the possibility to establish the presence of flat directions of the black hole 
potential. Classically, the flat directions are not fraught with any problem. Namely, the 
criticality condition of Vbh might not flx uniquely the values of the scalar flelds, and some 
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combinations of the fields might remain free, while the value of the entropy does not depend 
on these combinations. The problem might be hidden in quantum correction terms which, in 
general, might destabilize the flat direction (transforming it into a saddle critical point). 

The presence of the additional terms - whether of quantum or classical origin - in the 
action ([T]) leads to a modification of the black hole potential; there arises the so-called effective 
black hole Veff potential, whose extreme value is equal to the entropy of the black hole and 
whose critical points give the values of the scalar fields at the black hole horizon. For the case 
of model ([1]) the effective potential K// is equal to the black hole one Vbh- 

An alternative approach - the entropy function formalism - for calculating the entropy of 
a black hole with higher order derivative corrections was proposed in [6l [12] . It may also be 
applied to calculate the entropy of a large class of systems with additional degrees of freedom, 
e.g. rotating black holes, black strings, multi-center black hole configurations. Since we are 
also interested in stability issues, we first relate this approach to the black hole potential 
one and derive the concept of the Hessian matrix in the entropy function formalism. Let us 
recall that the latter is based on putting the theory on the near-horizon background geometry 
and constructing the so-called entropy function which can be easily written down once the 
Lagrangian is known 



2tt 



where stands for all fields but the moduli. In this formalism the black hole entropy is given 
by the value of the entropy function at the horizon. The horizon values and 0^ of all fields 
are determined by resolving all criticality conditions simultaneously 
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Since the entropy function formalism is more general, it encompasses the black hole potential 
approach. One can obtain the black hole potential from the entropy function as follows: 



VBHi<P)=n''£iE\<P),<P) 



(6) 



where E\(l)) denotes the solution to the equations dS/dE^ = in terms of the moduli. Of 
course, one can find the black hole potential, provided it is possible to resolve this equation 
explicitly. For example, in order to reproduce the black hole potential ([2D, one should fix the 
near horizon geometry as AdS2 x S"^ 



ds'^ = -vi r'^dt^ y] + ^2 [de"^ + sin^ 



(7) 



resolve Bianchi identities for the electromagnetic fields 
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-e'^(r). 



^23 =P sin 6* 



and substitute the explicit dependence of E^ = {e^,vi,V2) on 0° 
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back into the entropy function. 
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In many cases one cannot resolve the criticality conditions (jS^) and hence derive the black 
hole potential, therefore the entropy function formalism is more general. In order to make it 
self-contained one should be able to answer the question about the stability, that is to calculate 
the Hessian matrix in terms of the entropy function only. It becomes important when one 
cannot deduce the exphcit dependence of on 0°. 

From the formula (jS]) it follows that the Hessian matrix is equal to 
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Supposing that one does not have the explicit form of E\(j)), let us try to express dE^ /d(f)°' 
through the entropy function. For this purpose it is sufficient to differentiate dS/dE^ = with 
respect to 0", considering E function of 6"". 



d'^S dE' 
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which yields the following expression for the Hessian matrix: 
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Notice, that this expression is quite simple to deal with. It requires the knowledge of the 
derivatives of the entropy function which are easily calculable. Amusingly enough, for studying 
the stability of the solutions, one does not have to be able to resolve equations d£/dE^ = 
in terms of the moduli in order to reconstruct the black hole potential. This issue becomes 
important when considering higher order derivative corrections. 



3 Higher order corrections 



Now let us expand the considerations we made above by adding higher order derivative terms 
to the action ([1]): 
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The last term corresponds to higher order derivative corrections coming, for example, from 
the a '-expansion of the heterotic string action [13] ; its form will be specified later. 
The entropy function corresponding to the model ([U) acquires an additional term 
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which is defined as an integral calculated on AdS2 x S"^ background (171) with resolved Bianchi 
identities (IHD 



SH{vi,V2,p,e,4>) = -7 / dOdipJ - g Ch- 
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Since the contribution from the term 8h is supposed to come from the perturbative expansion 
of the superstring action, we present it in the following form: 

8H = a'8i + a'^82 + 0{a'^). (13) 

Our goal is to construct the black hole potential corresponding to the action ([U). In what 
follows we reserve the term "black hole potential" for a potential corresponding to the classical 
action ([1]), and the term "effective black hole potential" - to the action ([U) with higher order 
derivative corrections. 

As it was demonstrated in the previous section, to obtain the effective potential one should 
eliminate in the entropy function all fields except the moduli, using their equations of motion. 
For this purpose we represent the entropy function in the form 

8{e, V, 0) = 8o{e, 0) + a '^i(e, t;, 0) + a '^82(6, v, 0) + 0(a (14) 

and for simplicity introduce again a notation = {vi,V2,e^) for the "superfluous" fields. It 
can be easily shown that the effective potential, up to the second order in a ', is then given by 



K//(0) 



1 

TT 



8{E',r) 



a 



/2 



iJ) 



_i d8i d8i 
dE' dE ' 



+ 0{a 



H 



IJ 



J E=Eo 



dE'dE-' 



(15) 



where Eq stands for the "classical" solutions ([9]) for the fields f 1^2 and in terms of the moduli. 
The matrix has a following block form: 
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with Vbh and given in eqs. ([2]) and (Q. Therefore we may represent the effective potential 
in the form 

''%2(i^o(0),</>) 

(17) 
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with the right hand side, obviously, being calculated on the classical solution ([9]). 



3.1 Explicit form of the corrections 

Now we try to make our considerations more specific by defining the form of the corrections. 
Although we are interested only up to a'^ order corrections to the effective potential, the 
formulae of this section might be easily written down for the corrections of any order. 

To the a order, the possible corrections allowed by dimensinal analysis have the following 
schematic structure: 

where the Lorentz indices are supposed to be contracted with a proper number of inverse 
metric g^"^ in all possible ways. This expression induces the following form of the first order 
correction to the entropy function: 

£1 = ^B? + £rf'^ + ^F*- (18) 



6 



Calculated on the AdS2 x background ([7]) and with the field strength given by eqs. 
terms composing Si are then given by 
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where all a*-"^ are functions of 0°. 

Very schematically we present, as well, the second order corrections to the entropy function 
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Here the coefficients are assumed to depend on the scalar fields (j)"-, but their dependence is not 
written explicitly for the sake of simplicity. The function S2 f|T51) . in turn, is given by the sum 

^2 = ^_R3 + ^_R2F2 + ^_RF4 + £^F6- 

Each of the functions of eqs. (fT9l) and (!20|) turns out to be an eigenfunction of the operator 
appearing in the expression for the effective potential ([1] 
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This allows us to rewrite the effective potential fll7p in the following form: 
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where, as before, the right hand side is supposed to be taken on the solution (l9l). 
3.2 Corrections to the value of the entropy 

In the previous section we derived the effective black hole potential which encodes the informa- 
tion about the entropy of a black hole. Now we are going to calculate the value of the entropy 
up to the second order in a ' . So, we should extremize the effective black hole potential with 
respect to the scalar fields 

dVeff 



0. 



(21) 



In order to resolve these equations perturbatively let us expand the effective potential and the 
scalar fields in series on a ' up to the second order 



V,ff = Vo + a'Vi + a "V2 + 0{a '^), 



'20^ + O(a'=^). 



(22) 



7 



The expansion of Veff is nothing but a concise form of the eq. fll7l ). Being written in full details 
it gives 
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where all £^-terms are calculated on the "classical" solution (Q. After substitution of the 
expansions (122!) into the criticality condition (!2T|) one immediately derives 
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where the subscript zero means that the corresponding expression should be taken upon 0" = 
(pQ. Although, in order to obtain an extreme value of a function up to a order, it is enough 
to know a solution to criticality condition up to a order, we wrote down as well an 
equation fl24b ) defining the second order solution It relates to a subtle effect that we 
illustrate on a simpler example. 

Let us suppose that we are interested in calculating the value of the entropy up to a 
order. In this case we need to know only a "classical" solution 0g to the eq. (12^). so that 
one might think that eqs. fl24b ) are of no importance, since they define the value of the first 
order solution The subtlety comes out when the "classical" black hole potential Vq has 
a flat direction. In this case the matrix of the second derivatives (which is nothing but a 
Hessian matrix) becomes degenerate and not all of 0q might be determined from fl2^ ). The 
undetermined scalar fields appear then in the right hand side of fl24b ). This fact might make the 
system of eqs. fl24b ) to become inconsistent, since the Hessian matrix is degenerate. Considering 
this possibility as quite unphysical, we impose the condition that fl24b ) be consistent. This 
condition might fix some of the previously undetermined If the eqs. fl24b) turns out to be 
consistent identically, it means that the symmetry of the first order corrections to the effective 
potential coincides with the symmetry of the "classical" effective potential. 

When one is interested in corrections up to a'^, then one has to check also that fl2^ ) is 
consistent. The extreme value of the effective potential (that is, the entropy in fact) is then 
given by the expression 



eff.extr 



Vo{(j)o) +a'Vi{(f)o) +a 



12 



(25) 



Let us analyse this expression. First of all, the black hole potential ([2]) is a homogeneous function 
of degree two in the charges and q^. Then, the "classical" solution 0o is homogeneous of zero 
degree on the charges. It immediately follows from the fact that Vbh is homogeneous and 0o 
is a solution to a homogeneous equation (jl]). In order to write down the dependence of the 
entropy on the charges, let us for simplicity combine them to form a symplectic charge vector 
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Calculating the expressions Si^2 on the "classical" solutions according to formulae ffT^ and fl2U]) 
and taking into account eqs. (1231) . and (l25l) . one gets 
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The coefficients A, A^^, A^^-^^ and 5, B^^, . . . are composed of the previously introduced func- 
tions and u^^ and in general depend on the "classical" values of the scalar fields that 
are homogeneous of the zero degree on the charges. The explicit form of these coefficients 
depends on the model one considers and, if one considers higher dimensional theory, on the 
possible compactifications to four dimensions. Nevertheless, the structure of the contributions 
coming from the higher derivatives remains the same. 

The generalization of this formula for higher order derivative corrections is straightforward. 
We see that really all the series is built out of two quantities: the classical value of the entropy 
(proportional to Vbh), which is quadratic in the charges for the extremal black holes in D = 4, 
and the classical values of the scalar fields, which are homogeneous of degree zero. 

Formula (!26|) is an agreement with the result obtained for the one-loop correction to the 
Bekenstein-Hawking entropy [H] when only the term is present in the action ([TJ). 

As we have noticed before, the coefficients A and B (with any number of indices) in fl26p 
depend on the classical values of the scalar fields. If the black hole potential Vbh has a 
fiat direction, this means that not all of the scalar fields get fixed. This is not a problem 
for the Vbh, since its value does not depend on the scalars which are not fixed. But the 
first order correction might depend on these not fixed values. This means that black hole 
entropy depends on the values of the non-fixed scalars which can change continuously. This 
fact threatens the attractor mechanism paradigm and can induce a violation of the second 
law of black hole thermodynamics. Really, such a violation does not occur and the attractor 
mechanism paradigm is preserved. As it has been already explained, once a fiat direction is 
present, the system of eqs. (12^) becomes degenerate and to possess a solution its right hand 
side should satisfy some consistency condition, which might fix some of the previously free 
scalar fields 0q. If not all of the scalars 0q get fixed, then the group symmetry of the higher 
derivative corrections is a subgroup of the symmetry of the black hole potential, hence the 
higher order corrections will not depend on these non-fixed scalars. 

Looking at the eq. fl26|) one sees that higher order corrections become singular when con- 
sidering small black holes. Once Vbh is equal to zero, the perturbative expansion is not valid 
anymore. And there arises a question whether quantum effects might cure small black holes. 
This will the subject of a forthcoming analysis, whereas for the time being we limit ourselves to 
just some remarks. In realistic models the singularity might be removed when along with 5*0 = 
the coefficients A and B tend to zero as (p"' = 0g in a way such that the corresponding fractions 
in (l26l) remain finite. 



4 Application to the stu model 

In order to shed light on specific features of the effect of quantum corrections in the case 
of small black holes or when flat directions are present, we are going to revisit some well 
known examples where both such features are present. In particular we consider the stu model 
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with higher order corrections stemming from compactification of the heterotic string theory 
to four dimensions [71 [10]. The classical non-BPS solution for this model possesses two flat 
directions [15] in the axionic sector of its moduli space. In addition, this example allows us to 
investigate small black holes. 

In order to derive the form of the higher order corrections to this model we take the action 
of the heterotic string model [13], with all a' order terms calculated, and truncate it to six 
dimensions. Being afterwards reduced to four dimensions it corresponds to the stu model with 
higher order derivative terms. 

We perform our computations directly in six dimensions since dimensional reduction of 
the a ' corrections to four dimensions might be a topic of independent research. Moreover, 
working in six dimensions makes it easier to take into consideration the gravitational Chern- 
Simons term [71 [10] . 

We will follow the strategy given in [10] and, in order to simplify the comparison of the 
results, we will mostly follow the notations by [10]. The difference of our consideration from 
that of [TU] is that we preserve all the scalar fields - axions and dilatons - present in the theory, 
since the axion fields turn out to play an important role when it comes to the question of 
stability. 

The action [ISj we start with is the bosonic part of the heterotic string theory action 
truncated to six dimensions. It describes the coupling of a two-form field B^^ and a dilaton 
field $ to six-dimensional Einstein gravity 
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where 
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Here the term AC^ is of the forth order in the curvature R^^^^^. The field strength H^^^ 
includes also a gravitational Chern-Simons term 
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^PQR = J^P (Kn^'^mr) + I H^'pn^qH^'mr - \ R^'^'pgH' 



rRMN 

^ - P \' QJV — MR/ ' g — PN ■ V — MR ^ ' 

+ — p^H"" ^^^H'' + antisym. on P, Q, R. 



(29) 
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As prescribed in [10], one may dualize the field strength H^^^ into a new one K^^^^ which is 
an exact three-form 

^PQR = '^P^QR + ^Q^RP + ^h'^PQ' 

obeying, obviously, the Bianchi identities. To this end, one adds the term 



into the original action (1271) and eliminates the field strength Hp^p. We postpone the elimination 
of H„^„ for a while, for a reason to be clarified later. 
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Dealing with the Chern-Simons term is a httle bit tricky [TU] and consists in maintaining the 
Lorentz covariance of the Lagrangian. We succeeded in singhng out a manifestly covariant part 
of the Lagrangian by adding total derivative terms, in such a way that the second derivatives 
of the four dimensional electromagnetic potentials originating in the six dimensional metric 
disappear. 

In order to be able to interpret the model (127|) as an ancestor of the stu model, one should 
clarify the field content of the theory. First of all, the six dimensional metric tensoio G^,^^ when 
reduced to four dimensions produces a four dimensional metric tensor gf^jy, two vector fields Ajji"^ 
and three scalar fields - two dilatons ui, U2 and one axion c. Then, the two-form potential C^^j^ 
produces a four dimensional two-form C^^J , two vector fields Aj^'^ and one axion Cmn = bemn- 
In four dimensions a two-form C^J is dual to a scalar; this duality gives rise to another axion a. 
Finally, the original scalar field $ gives rise to a four dimensional dilaton Ug- So, we end up 
with a four dimensional metric tensor, four vector fields and six real scalars - three axions a, b, c 
and three dilatons ui,U2,Us (more precisely, exponentials of the dilatons). 

The entropy function formalism in this case is based on putting the theory on a six- 
dimensional background 

rfsjg) = dsf^^ + Gmn {dx"^ + A^-'^dx^') (dx" + A^-^dx") , (30) 
where the two-dimensional metric tensor Gmn is parameterized as follows [16j: 



„ _( u\-\-C^u\ Cu\ \ 



— CU2 U2 



The connections ^ are chosen in the following form: 



Al={2re\ 0, 0, o) , = (2 r e^, 0,0,-^ cos 6). (32) 



Hence, in what follows they will give the so-called magnetic configuration corresponding to 
the DO — DA branes. The other two vector fields coming from C are equal to 



MN 



4^ 



f -re3,0,0,0) , = f -re^O,0,-^cos^ V (33) 

V g ' ' ' y ' ^ V 8 327r ) ^ ' 

The electric potentials e^- are dynamical fields; therefore we do not put them equal to zero ad 
hoc, as their values are to be fixed when minimizing the entropy function. 

In the entropy function formalism one should fix the values of the scalars at the horizon of 
the black hole. Thus we parameterize the dilaton field as follows jlO] : 

Us 

From a four dimensional perspective the field strength K^^^^^ is related to the vector potentials 
and axions in the following way: 



Us 



"'^we assume the splitting of the six dimensional indices M,N,K,... — 0,1,..., 5 into four dimen- 
sional fi,v, p, . . . = 0, 1, 2, 3 and two dimensional m, n, fc, . . . = 4, 5 ones 
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and, therefore, its non-zero components are given by the expressions 



Koi3 = -^p\me'-e^)cose, Kou = ^ + 2be'', Koi5 = ^~2be\ 

^023 = -T^ (166eV + eV)'^sin6', ^^234 = -— sin6', K235 = -^7^sin6'. 

As it has been mentioned above, the field strength -ff^,^^^ should be eliminated by means of 
its equations of motion. Since we perform our calculations up to a order, one has to resolve 
these equations of motion up to the a'^ order, what seems to be quite sophisticated due to 
the nonlinear nature of the action fl27|) . In order to eliminate if^^^^, we perform the following 
trick [7J: we use the Ansatz (dictated by the symmetry of the problem) for its non- vanishing 
components 

Hoi3 = hicos6, iifoM = ^2, Hqi5 = -hi, 

(35) 

-f^023 = h^r sin 9, H234 = sin 9, H235 = h sin 9 

and include the fields hi in the set of dynamical fields, with respect to which the entropy 
function is to be minimised 



d£ _ dS _ d£ _ d£ _ 
dhi~ W ~ d^i~ d^2~ 

da db dc dui du2 dug 
The entropy function acquires the form 

8 = 2tx [e^qi + e^gg + - - IQbhiix - 326/i47re^ - ?,2bh^Tie^ - 2h^'Ke^ 

, ^, 4aeV h2P^ , aeV ^3^5 , h3hie^Us hl{e'fus , /la/isc^u^ 
+z/i47re — 1 — 1 h 



(36) 



47r 2 47r lQv2 ^V2 4^2 4^2 

hJi^e^e^Us h\{e^Yus ^ h\usVi ^ ch/^h^UgVi ^ h^UgVi ^ c^h^UsVi 

2t>2 4t>2 IQU1V2 8M|f2 16^2^2 lQu\v2 

u\{p^)'^UsVi h\usV2 u\e\usV2 (?u\e\usV2 cu\eie2UsV2 
64:Tc'^V2 IGufvi 4t>i 4t>i 2vi 

u\e\UsV2 (?Tl'^h\UsV2 Ti'^h\UsV2 CTlhih2UsV2' 



(37) 



4t>i Auiip^Yvi Aulip^Yvi Auip'^vi 

Here, for the sake of brevity, we avoid reporting the explicit expression we derived for the 
first a' order corrections. When minimizing this expression, we perform an expansion of all 
dynamical fields over a', i.e. 



a 



a(o) + a 'a(i) + . . . , b = 6(0) + a 'b(i) + . . . , etc. 
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In this way one gets "classical" solutions to fl5B]) 



,(0) 



(0) 



,(0) 



Ho) 



C(0) 




(3^ 



\ e"3(l + e2 

Here the three real parameters at are restricted by one condition [17] 

«! + 02 + as = 0, 

so that the presence of two independent unconstrained parameters indicates the presence of 
two flat directions. 

The charges above are connected with standard stu black hole charges by the following 
t r ansf ormat ions : 



gi = 4 7rQo, 



p 



p3 



p 



One can check this by making a transformation from 5*0(2, 2) to {SU{1, l)/f/(l))^ basis pTlfTS] . 

Naive substitution of the zero order solution into the entropy function makes it dependent 
on these parameters. This effect appears to contradict the attractor mechanism paradigm. 
However, a subtler analysis shows that these parameters get fixed and the correctness of the 
attractor mechanism is restored on the quantum level too. 

To illustrate how the dependence on these parameters drops out, we present in details the 
first order correction to the entropy. So, the "classical" solution fl38l) yields the value of the 
entropy up to the a'^ order 



3 + 2e2°2 + 3e4"2 
(1 + e2°2)3 



p'^qi 



+ 0{a 



i2\ 



Let us note that even including axions the non-BPS solution the first order correction to the 
entropy does not contain any contribution form the Chern-Simons term, in an agreement 
with [10]. Sticking to the first order approximation for the entropy, the only thing to know 
about the first order solution to the eq. fl36l) is that it exists. In general, this is not guaranteed 
automatically, since the matrix of the second derivatives of the entropy function is degenerate. 
In the case under consideration the first order solution exists if 



02 = 0. 
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In this case the first order solution acquires the following form 

32e2°i7rgi 



(1) _ aip"^ y - P'^P^QiQs 



(1) 



h 



6(1) 



C(l) 



u 



(1) 



(l + e2 



^1 + e^^i) ai(p2) - 96 (-1 + e^^^) tt^ 



, ^5 



(1) 



-1 + e2"i)V 



(l + e2°i)>2 



16(l + e2"i)7r3.. 
;i + e2"i) ai(p2) Vgs + 64 (-1 + e^^i) vr^/ - pV?i?3 



47rp^ 



8(l + e2"i)7rVp^?3 
3 47rp^gi 



p^p^qiq3 



93 \/ -PV9193 



^(1) 



0, .«=.«=0, 



, , -128e2"i (-1 + e^"0 nbip^sJ -p^P%q3+ 

8 (1 + e2"i)3 vrp2(p4)2yr|rg3 L V ; y3v y ^3 

(^1286^"^ (l + Ge^"^ + e^"^) vr^gi + (-1 + e^"^) (l + e^"^)^ 
:i + e2"i) 6i(p2)2^2 _ (_i + e2°i)7rV-pVgig3 



-1 +e 



(1 + e2"i)2 (p2)%/fV-pVgig3 



16 ri + e2"i)V (-p2p4gig3) 3/2 

( -1024e2''i (-1 + e^^i) Tr^gi + 4 (l + e^^i 
1287r^A/-p2p4q,^gg 



1024e'"^ (l + e2"i j -nh^q^^ -p^p^q^q3+ 

3 



'P'^p'^qiqs 



{p^fql 

Let us note that the two scalar fields (in our case these are a(i) and remain undefined in 
this approximation. With vanishing a2 the first order correction to the entropy does not depend 
on any free parameter anymore and the second order correction acquires a much simpler form, 
so that we may write it down as 
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2„ / 



p p qiqs + I6n a 



-—^ 167r a — , - + 0(a^) 

P^3 P^qz J - p'^p^qiqs (1 + e2"i f 



Performing analogous steps to second order, one can check that requiring the existence of 
the second order corrections to the solution of eqs. (136!) yields ai = 0. The solution we obtained 
in this way, corresponding to = 0, coincides with the axion free solution given in [71 [TO] . 
The only difference is that we derived this solution keeping the full axion dynamics, without 
truncation. It is only on the horizon that the axion contribution vanishes. 

The value of the entropy up to the second order is then given by 



S 



2„ / 



p p qiq3 + 167r a 



pSi 190 4/2 Ai 
- l/87r a 



p'^q3 



p'^qs 



p'^p^qm 



+ 0(a 



i3\ 



(39) 



One can see that this expression shares the general features pointed out in ([26 
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Exact solutions 

Despite the cumbersome and entangled structure of the above non-BPS perturbative solutions 
to the extremization condition of the entropy function, one can cast them in exact form|§ 



Vi = V2 

h2 = - 



47r2 ' 



2^/ ' 



V-p^giP^gs /, , 327r 



2a' 



u ^3 

3 _ yJ-p^qiP'^q^ 



(40) 



47rgi 



Ml 



9i 



r ,/i + 



'93 
^2 ' 



Us 



4vrg3 



1 I 327r2o' ' 

-L T 



P 93 

1 



for the charges 



qi < 0, g3 > 0, > 0, 0. 



(41) 



Once the close form of the solution is obtained, one can see that the perturbation expansion 
that we were performing is valid when 



32 7r2a' 



< 1. 



If one fixes the value of a then one may understand this formula as a condition on the charges 

b'gsl » 1, (42) 



when classical effects become dominant. As we will see later, the condition ( H2l) fails for small 
black holes. 

In the paper [7] another exact non-BPS solution was found, which in our notations rea dP 



Vi = V2 
h = - 



P% 

'4vr2 



2^ 



327r^a 

p'^q?. 



P^qip'^qs 



V-p^qip'^qs 

47rgi 



327r2a' 



p qs 



Ui 



U2 



2 



p 



\ 



327r2a' 

2 ' 

p q3 



q3 (. 32n^a'' 



p^qz 



- p^qip'^qs 
Anqs 

Snp^qi 



327r2a' 



P^qiP^qs I S2fa' 

V p^q^ 



(43) 



which is valid for the following charges: 

qi < 0, q3< 0, p^ > 0, < 0. 



(44) 



^only non- vanishing fields are written down 
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Both solutions correspond to the same value of the entropy 

S = 



1 + 



327r2a' 



(45) 



In the classical hmit a' = the Hessian matrix of the solutions (l40l) and (l43l) has the following 
eigenvalues: 



Ai = A2 = 0, A3 = 2|p^|W- 



As = 



p'^p^qi 



P^qi 



{prqi + 2hQ{p'fqi + 



(46) 



167r^ 



93 



When turning on the quantum effects, one finds that the eigenvalue Ai remains equal to zero, 
while A2 acquires a positive correction 



A2 = 327r"|p4k - 



p4 



(/)2 + AmQu'^ql 



p^q^qs 167r2(p4)2g2 ^ (^2)2(^4)2 _^ 40967r2g32(p2^2 



(47) 



so that one can say that the solutions are "stable" . 

To complete the list of solutions, we present as well two additional solutions. One of them, 
found in [7j, is a BPS solution (see footnote [2] on p JTSl) 

327r2a'\ 



Vi = V2 



P% 
47r2 



1 



p2g3 



327r^o 



2„/ 



2Tip^ 



1 - 



647r2c 



q3_ 

27r 



327r2a'' 

p'^qs 



p'^p'^qiqs 



647r2Q;' 

p q3 



(4J 



Ml 



?1 



327r2a' 



1 - 



U2 



4vrg3 

p2 ) 



647r2c 



327r2a' 

p^qz 



Stt 



647r2a' 



and is valid for the charges 
The entropy is then given by 



gi > 0, gs < 0, p^ > 0, / < 0. 



S 



1 + 



647r2Q;' 



\ Ip'^Qal 



(49) 



The other solution acquires no a ' corrections 



Vi = f 2 



P^g3 
47]-2 



1 p r93 



47r V p'^qi 



ho = - 



qi_ P Qs 
27rVpV 



q3_ 

27r 



U2 = 



p^qz 



(50) 



16 



It is valid for the charges 

qi > 0, gs > 0, > 0, / > 
and the corresponding entropy is equal 



S = ^p%p'^qs. 

Despite the fact that it has no quantum corrections, this solution is non-BPS but with vanishing 
central charge. These two solutions are stable in the classical limit, and we suppose that 
quantum effects do not spoil this feature. 



Small black holes 

Out of the solutions derived in closed form in previous section, only two admit small black 
hole limits, i.e. a 1/2-BPS f l48p and a non-BPS Z fH5]) . solutions. The remaining two 
solutions are non-BPS and do not contain small black holes (one of them, i.e. the non-BPS 
Z = solution, remains even unaffected by quantum corrections). By definition, a small black 
hole has vanishing classical entropy. In the solutions reported in the previous section, the small 
black hole limit corresponds to 

ga = 0. (51) 

When this limit is considered, only the two solutions fHS]) and fH51) remain regular, while for the 
remaining two the radii of the AdS2 and S'^ spaces go to zero. Furthermore, the non-BPS Z = 
solution (!50|) has vanishing entropy, whereas the non-BPS Z ^ solution (140|) has the same 
(non vanishing) value of the entropy as the solution (H3l) . 

One can easily understand that the perturbation theory over the parameter a ' fails for small 
black holes. Indeed the genuine parameter to make a perturbative expansion is for small 
black holes and l/gs for large ones. 

Taking the limit (ISTjl in the solutions P3l) one gets the non-BPS small black hole solution 



(52) 




a' ' 



with the corresponding entropy 



S = ATi^2a'p^qi. 

Calculating the Hessian matrix on the solution fl52l) . one obtains the following eigenvalues: 

(p2)4(p4)2 ^ 36(87r)6(p2)2 ^ Q{ATif{p^f {a'f 



Ai = 0, As 



512y27r3(p2)2 /p4^^^/ 



A3 = — —pJ2a'—, A4 = , , A5 = 0, A6 = 0. 

3 V Jp%a' 



(53) 



One sees that A24 < and A3 > 0, which means that the solution (l52l) is not stable. 
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Now, taking the limit flSlI) in the solutions fl48l) one gets the BPS small black hole solution 



Vl = V2 = 8a', h2 = 2, 



3 -p-gi 
with the corresponding entropy 




(54) 



In this case two of the eigenvalues become zero, three of them become positive and one becomes 
negative 



0,1 + yaf + bi ai - y'af + 6i 02 + y 03 + &2 ^2 - y ^2 + ^2 

Cl Cl C2 C2 

where 

ai = 25((/)y)2 + 1145(47r)6(256(p2)2 + {py)a'\ 02 = 7(/)V^i + 2887r\p''fa' - Sevr^a'^, 

61 = -1294(8^)^((/) V)^(256 (p^)^ + {py)a", 62 = 121607rV)^-(p2)y + gvr^a a', 



Cl = 2588(27r)'^(p^)Y-p''9ia', C2 = SSvry^ -pVo 

We can then conclude that, amazingly enough, for small black holes, all found solutions - BPS 
and non-BPS - are not stable. 

5 Conclusion 

In this paper we considered a correction to the black hole entropy due to the most general higher 
order derivative terms in the Einstein-Maxwell action. We demonstrated that the general form 
of the corrections to the entropy is in agreement with previously found results [H]. Provided 
that the perturbation expansion over classical solutions is valid, the form of the correction is 
completely determined by the classical value of the entropy 5*0 (a homogeneous function of 
second degree in the charges) and the classical values of the moduli (a homogeneous function 
of zero degree in the charges). 

The fact that the subleading corrections are singular in Sq drops us a hint that small 
black holes are purely quantum objects. In fact, the considered example of the stu black hole 
illustrates how small black hole solutions are singular in a', so that the standard perturbation 
expansion fails and one should find another small parameter to carry out perturbation theory {q^ 
rather than l/qs)- Notice, however, that the inadequacy of the parameter a' as a perturbative 
parameter for small black holes is valid in general, not only for the considered example of 
the stu black hole. The only feature that will depend on the model is the explicit form of 
the genuine perturbation parameter corresponding to the small black hole limit. It would be 
interesting to find a model independent definition of the parameter, generalizing its specific 
realization (!5T|) in the model treated in this paper. 

We argued that the attractor mechanism remains valid in the presence of higher order 
derivative corrections, despite the fact that in the presence of flat directions in the classical 
black hole potential, the corrections to the entropy apparently might depend on the undefined 
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moduli. This is made possible by requiring, to each perturbative order in a the fulfillment of 
a consistency condition, needed for the existence of a solution to a degenerate equation system, 
hence fixing some of the undefined moduli fields. If not all of the scalars 0o get fixed, then the 
group symmetry of the higher derivative corrections is a subgroup of the symmetry of the black 
hole potential, hence the higher order corrections will not depend on these non-fixed scalars. 

We established a relation between the black hole potential approach and the entropy function 
formalism. The black hole potential can be deduced from the entropy function by eliminating 
"superfiuous" fields from it. In the presence of higher derivative corrections this procedure 
might be quite difficult, if not impossible. Even when deriving the black hole potential proves 
difficult, the issue of stability of the solutions can be studied quite easily, entirely within the 
entropy function formalism. In this way the stability of the solutions was studied in this paper. 
We believe, however, that the validity of proposed method is quite general and its applicability 
goes beyond the considered model, and even beyond the considered class of systems - black 
holes in four dimensions. Hence, it would be interesting to further apply the method to study 
the stability also of extended objects, such as black strings and black branes, as well as rotating 
black objects, considered not only in four but also in higher dimensions. 

In this paper we chose a specific model, i.e. the N = 2 d = 4 supergravity model with stu 
prepotential, since it is known to possess fiat directions in the non-BPS branch fT5]. In order to 
derive a form of the corrections pertinent to this theory, we started from the action of heterotic 
string theory [T3j with all a' corrections included [7j. Consequent compactification down to 
four dimensions reproduces a stu model with all necessary corrections included. The most 
interesting features for us reside in this case in the sector of the axion fields; therefore we did 
not make any truncation over the axions. 

Apart from the known solutions of the obtained model, we found two more non-BPS ones. 
It is noteworthy that one of these solutions acquires no a' corrections and it is a non-BPS 
solution with zero central charge. 

We investigated the small black hole limit of the above mentioned solutions and found that 
the solutions found by us are singular in such a limit, while the previously known solutions 
become unstable in the sense that the corresponding Hessian matrix acquires positive, negative 
and zero eigenvalues. It is of interest to investigate whether this property of small black hole 
solutions is general or just a peculiarity in the considered stu model. 
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